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The approach in a theory of collective excitations in hot nuclei exploring the formalism 
of thermo field dynamics and the model Hamiltonian consisting of a mean field, the BCS 
paring interaction and long-range particle-hole effective forces is reexamined. In contrast 
with earlier studies it is found that a wave function of a thermal phonon is depended 
not only on the Fermi-Dirac thermal occupation numbers of Bogoliubov quasiparticles 
consisting the phonon but on the Bose thermal occupation numbers of the phonon as 
well. This strongly affects a thermal phonon couplings due to renormalizing of a phonon- 
phonon interaction and enlarging the number of thermal two-phonon configurations 
coupled with one-phonon ones. Moreover, it is shown that the formulation of the double 
tilde conjugation rule for fermions proposed by I. Ojima is more appropriate in the 
context of the present study than the original one by H. Umezawa and coworkers. 



1. Introduction 

The present paper addresses some aspects of theoretical treatment of collective 
excitations of a hot nucleus within the thermo field dynamics (TFD). 

The TFD 1 z is known as a powerful tool in studying many-body problems at 
finite temperatures. There exist numero us TFD applications in condensed matter 
and high energy physics (see e.g.^EEEJ). In spite of several attractive properties 
the TFD formalism seems to be less popular in the nuclear theorist community. 
During the last 20 years only a dozen works dea ling with the TFD applicatio n to 
nuclear structure problems were publishedP ^ | 8 | 9 | 11 ' 1() | 12 1 13 1 14 1 15 1 16 1 ^ 1 18 1 19 1 

An important step in application of the TFD formalism for extending the stan- 
dard nuclear theory methods (the HFB approximation, RPA, a boson exp ansion 
technique) to finite temperatures was made by T. HatsudaP injBIIIIISIISEI differ- 
ent approximations going beyond the thermal RPA were formulated. O. Civitarese 
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and coauthorJ^ES applied the TFD to problems of treating nuclear pairing in the 
BCS and the random phase approximations as well as quasi-continuum (or reso- 
nant) single-particle states at finite temperatures. Moreover, the TFD was used to 
prove the Bloch-Messiah theorem at finite temperat uresP and construct number- 
projection methods for the BCS pairing in hot nuclei PSliHI 

The methods and approximations developed in the papers cited above were not 
applied to numerical calculations of any nuclear properties. Only simple so lvable 
models were used to demonstrate the new achievements. At the same time, irfffilE] 
the TFD was combined with the Quasiparticlc-Phonon nuclear model 

(QPMjZQl 

to formulate an approach allowing one to treat microscopically a coupling of dif- 
ferent excitation modes in hot nuclei. This coupling is known to be responsible 
for the spreading width of giant resonances in cold spherical nuclei, and the TFD- 
QPM approach was used to analyze a thermal behavior of a giant dipole resonance 
widthM^ 

It appeared that in contrast with earlier studies of the proble (see alsc 
based on the Matsubara Green function technique and the Nuclear Field Theory, 1 ^ 
the quasiparticle-phonon interaction at finite temperature in the TFD-QPM ap- 
proach did not depend on the thermal occupation numbers of phonons. Only the 
Fermi-Dirac thermal occupation number s of no ninteracting BCS quasiparticles ap- 
peared in c orre spondin g fo rmulae of Refs! ^|18 | g ome aspects of this difference were 



discussed irffSl ( se e alscES). 

Recently, studying weak transitions of the Gamow- Teller type in hot nuclei^ 
within the TFD we found some inconsistencies in the TFD-QPM approach off^HSl 
related to structures of thermal vibrational phonons. Moreover, new points in the 
TFD-QPM app roach a ppear if one adopts changes in the TFD formulation proposed 
by Izumi Ojima^EZl that concern the alternative formulation of the double tilde 
conjugation rule. 

The re asons given above compel us to reexamine the TFD-QPM approach 
O j-10 | 12 | 1 8| t a quasiparticle-phonon coupling at finite temperatures. This is the 
main aim of the present paper. 

The paper is organized as follows. A brief summary of the TFD formalism is 
given in Sec. O This seems to be necessary because we u se here a variant of the the- 
ory which somewhat differs from the standard oneP^^In Sec. 02 the TFD formal- 
ism is applied to the nuclear Hamiltonian consisting of a mean field, the BCS pairing 
interaction and separable particle-hole forces, i.e., to the QPM HamiltonianPSl In 
Subsec. 13. 2[ the th ermal BCS pairing is considered. This part of the paper is close 
to that ofEHnHH a discussion of 7-transitions between a thermal vacuum state 
and thermal quasiparticle excitations is the only addition. In Subsec. 13.31 the equa- 
tions of the thermal random phase approximation (TRPA) are derived. It is shown 
why and how the thermal Bose-Einstein occupation factors should appear in the 
expressions for thermal q uasipartic le amplitudes of a thermal phonon wave func- 



quasipartic 
in 10 12 18 



tion. The points missed ir l u | 1 1 are analyzed and discussed. The new formulae 
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of the quasiparticle-phonon coupling term at finite temperature are evaluated in 
Subsec. 13.41 The summary and conclusions are given in Sec[U 

2. The TFD formalism 

We consider a system of Fermi particles at finite temperature and treat its sta- 
tistical properties within the grand canonical ensemble. In the standard statistical 
mechanics a heated system is described by a mixed state density matrix p which is 
the solution to the Liouville - von Neumann equation 

i% = [H,P\- (1) 

Here H is the Hamiltonian of the system under consideration with eigenstates |n) 
and eigenvalues E n (the chemical potential is included in H). The thermal average 
of an arbitrary operator A is given by 

((A)) = Tr[pA]/Tr[p] = £ e - E ^ T (n\ A\n) , (2) 

n n 

where T is the temperature in units of energy. The main idea behind TFD is 
to define a special state which is named a thermal vacuum |0(T)), such that the 
thermal average of A equals the expectation value of A with respect to this state 

((A)) = (0(T)\A\0(T)) . (3) 

To construct |0(T)), one should double the Hilbert space of the system by adding 
the so-called tilde states \n)U^ These tilde-states are the eigenstates of the tilde 
Hamiltonian H with the same eigenvalues as H, i.e. H\n) = E n \n). Thus, the 
Hilbert space of a heated system is twice as large as of the corresponding initial 
(cold) one. 

In the enlarged space covered by the states |n) £g) |m) two types of operators 
are acting - the ordinary say, Ai and the tilde ones Ai. There is a one-to-one 
correspondence between these two sets Ai <-> Ai. However, ordinary operators can 
change only ordinary states \n) whereas tilde operators change only tilde-states \n). 
There exist the following rules of tilde conjugation operation^ 

(AiA 2 ) = AiA 2 , 
( Cl Ai + c 2 A 2 ) = c* x A x + c*A 2 , (4) 

where ci, c 2 are c-numbers. The asterisk denotes the complex conjugation. The 
tilde operation is supposed to commute with the Hermitian conjugation 

(l)t = At , (5) 

Moreover, it is required that ordinary- and tilde-operators should commute or an- 
ticommute with each other 



[A 1: A 2 ] T = 



(0) 



September 8, 2008 8:59 WSPC/INSTRUCTION FILE dzhioev 



4 Authors' Names 



depending on their bosonic or fermionic nature. 

Introduction of the doubled Hilbert space of a heated system enables one to 
construct the above thermal vacuum |0(T)) in the following way: 



|0(T))=^e- B "/ 2T e M "|n)®|n)/ /^e-^/T, (a n £ K). (7) 




Hereafter, the thermal vacuum ([7]) will be referred to as the exact thermal vacuum. 
The thermal vacuum is invariant under the til de o peration. 

Till now our discussion follows the line ofM^ The new point appears in the 
definition of the so-called double tilde conjugation rule (DTCR). Originallyf^DTCR 
was introduced in the form 

A = p A A, (8) 

where pa — 1 if A is a bosonic operator and pa = —1 if A is a fermionic one. This 
DTCR form was used in all of the TFD ap plications to nuclear structure problems, 
i.e., i n . | 6 | 7 | 8 | 9 | 11 | 10 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | I 9| However, quite long ago the other form 
of DTCR was proposed by I. Ojima, u namely, 

A = A. (9) 

That is, the Ojima's version of DTCR does not distinguish between bosonic and 
fermionic operators. This form is based on the equivalence between the algebraic 
structure of TFD and that of the axiomatic statistical mechanics (c*-algebra ap- 
proach) established by I. Ojima. In the present paper, we will compare consequences 
of using the two DTCR versions. For that we define DTCR for a fermionic A in the 
following way: 

A = -a 2 A, (10) 

where a is equal to either 1 (A = — A) or i (A ~ A). Thus, the DTCR definition 
(flU)) includes both the above DTCR variants. 

The Schrodinger equation for a hot system in the doubled Hilbert space reads 

i^m,T))=nm,T)), (ii) 

since 7i = H — H is the time-translation operator and, therefore, the Hamiltonian 
of a hot system. The operator Ti is named the thermal Hamiltonian.^ The exact 
thermal vacuum ([7]) is the eigenstate of H corresponding to the zero eigenvalue. 
Equation (|lip describes the time evolution of a system at finite temperature and, 
thus, is the analog of the Liouville - von Neumann equation ([!]). 

The properties of elementary excitations of the system at T ^ are determined 
by H. Excitation energies of various modes at T ^ are the eigenvalues of H but 
not of the original Hamiltonian H and in general they depend on T. Moreover, 
any eigenstate of Ti. with positive energy has its counterpart - the tilde-conjugate 
eigenstate with negative energy. Following H. Umezawa we consider creation of a 
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tilde-state with negative energy as annihilation of a thermally excited state. This is 
a way to treat excitation and de-excitation processes in a heated quantum system 
within TFD. 

Whereas the dynamical development of the system is carried by the thermal 
Hamiltonian, its thermal behaviour is controlled by the thermal vacuum. The grand 
canonical value of an observable corresponding to operator A should be calculated 
as (0(T)|A|0(T)), i.e., after diagonalizing H. 

Obviously, in most cases one cannot diagonalize H. exactly and, thus, find the 
exact thermal vacuum and other eigenstates. Usually, one resorts to some approx- 
imations and then finds an approximate thermal vacuum state, e.g., the thermal 
vacua corresponding to the HFB or the random phase approximations. In the case 
that there appear several solutions in the given approximation, one should find the 
minimum of thcrmodynamical potential £1 to see which of them is realized. In TFD, 
is given aJ^2l 

Cl=(^ {T)\H-KT\^ (T)), (12) 

where ^o(T) is the approximate thermal vacuum and K is the entropy operator of 
the system. 

3. The Finite Temperature QPM 
3.1. The thermal QPM Hamiltonian 

In what follows we will use the microscopic Hamiltonian of the Quasiparticle- 
Phonon nuclear modePSl 

-ffQPM = H sv + i/p a i r + -ffph- (13) 

The Hamiltonian (|13p includes average fields of protons and neutrons 

4p-EE T ( E rA T )iv ( 14 ) 

r jm 

pairing interactions of the BCS type 

#pair = - J E Gt E 4imi4lmr a Hm5^j2m 2 ( a WT = (-l) J ~ m aj-m) , (15) 
t jimi 

.72 "12 

and effective multipole-multipole isoscalar and isovector forces 

ff Ph = ~\ E E (4 A) + P^ X) )Ml(r)M^{pr). (16) 

A/J. Tp— ±1 

The single-particle operator (r) reads 

M U T ) = ^ T {hrni\i x Rx(r)Y x ^9,cl>)\hm2)a] imi a j2m2 . (17) 

j2m 2 
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Here, an operator aj m (a jm ) is the creation (annihilation) operator of nucleon in 
a single-particle subshell with quantum numbers nljm = jm and energy Ej ; index 
t = n, p is the isotopic one, changing the sign of r means the interchange n <-> p, the 
notation ^ r implies a summation over neutron (r = n) or proton (r = p) single- 
particle states only. The parameter G T is the proton-proton or neutron-neutron 
pairing interaction constant; R\(r) is the radial form factor of the A-pole separable 
interaction, Y\ fl (9, cf>) is the corresponding spherical harmonic, Kq X ^ and k[ are the 
coupling constants of isoscalar and isovector multipole-multipolc interactions of A 
multipolarity. The value A r is the neutron or proton chemical potential (the Fermi 
level). 

To find an excitation spectrum of a hot nucleus governed by the QPM Hamil- 
tonian (fT3|) . at the beginning we should construct the thermal QPM Hamiltonian 
Wqpm 

Hqpm = -ffQPM — Hqpm, (18) 

where -ffqjpM = H sp + H pa i r + 7J p h is the tilde counterpart of TJqpm created by 
the tilde conjugation rules (j4]). Then we should diagonalize Hqpm- This procedure 
is quite similar to that used in the standard QPMpS i.e., at T = 0. The main 
difference lies in the doubled Hilbert space of a heated nucleus. 



3.2. Thermal quasiparticles 

The first step is diagonalization of part of the full thermal Hamiltonian, namely 
the sum of the two first terms Ti. sp + 7i pa ir which in the following will be referred 
to as the thermal BCS Hamiltonian Hbcs- To this aim, we make the Bogoliubov 
u, v- transformation from nucleon operators a\a to quasiparticle operators aJ , a 

a ]m ~ U j a ]m — v j a Jm , 

a jm = U j a jm - V i°7m » + V ) = !) ' ( 19 ) 

The same transformation with the same u, v coefficients has to be applied to nu- 



cleonic tilde operators a\ , a ?m , thus producing the tilde quasiparticle operators 



aj m and a 



jm jm 

Thermal effects appear after the second (or thermal) Bogoliubov transformation 
which mixes ordinary and tilde quasiparticle operators and, thus, produces the 
operators of so-called thermal quasiparticles 0j m ,Pj m and their tilde counterparts 

Pjm = X 3 a )m ~ Vy 3 ajm , 

P]m = X ^)m + m a 5m , + y] = 1) • (20) 

Note that in contrast with Refs!^- ^ | 10 | 12 | anc j man y others we include the factor a 
from the definition of DTCR (fTU|) to the thermal Bogoliubov transformation (|2"0)) . 

To find the coefficients u, v, we express the thermal Hamiltonian in terms of 
thermal quasiparticle operators (|20p and then require that the one-body part of the 
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thermal BCS Hamiltonian has to be diagonal in terms of thermal quasiparticles. 
This leads to the following equations for the u, v coefficients: 

-HO- 5 ^) (21) 

where Sj — -J (Ej — A T ) 2 + The gap parameter A r and the chemical potential 
A r are the solutions of the equations 

A - = ^E t ( 2 j + ^ - y>w> N * = E T ^ + l )^ x 1 + u M)> ( 22 ) 

3 3 

where N T is the number of neutrons or protons in a nucleus. Actually, the equa- 
tion for N T is not a consequence of the procedure described above but rather an 
additional demand that a number of particles of any kind in the heated system is 
conserved in average. 

With Uj , Vj from (f2~lj) the one-body part of the thermal BCS Hamiltonian reads 



Hbcs = H sp + H pail ~ Y^ e M^Jm - flrJjJ- (23) 



r jm 



Thus, 7Ybcs describes a system of noninteracting thermal quasiparticles and tilde- 
quasiparticles with energ ies Sj and —Sj, respectively. The vacuum for thermal quasi- 
particles is given b y I I 

\O(J3,0)) = exp{-^ / /2}e X p{ f 7*^^ 7 'at m St m }|0(a))|0(5)) , (24) 

r jm 

where |0(a)) and |0(5)) is the vacua for ordinary and tilde Bogoliubov quasiparti- 
cles, respectively. The operator Kf is the entropy operator. It reads 

&f = ~ E E i a ]m a jm ln Vj + a Jm a jm ln x ) } • ( 25 ) 
r jm 

We should stress that although the vacuum (HM)) is the eigenstate of the thermal 
BCS Hamiltonian (|23[) with zero eigenvalue it is not yet a thermal vacuum state in 
the sense of ([3]). To determine the thermal vacuum state corresponding to 7Ybcs we 
need to fix appropriately the coefficients Xj , yj . They can be found by minimizing 
the thermodynamic potential 

= E E r { £ ^ 2 + ^y) + x2 3 ln ^ 2 )}- (26) 

r jm 

Note that f2/ contains the ordinary operators H sp and H pa \ r but not the thermal 
ones 7i sp and 7i pa ir- As a result of variational procedure, we obtain 

-1/2 



l + cxp(!)p , X] = {l-y]) 1/2 . (27) 
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Thus, the coefficients yj are nothing else than the thermal occupation factors of 
the Fermi-Dirac statistics. 

The thermal quasiparticle vacuum ([24]) with the coefficients yj , Xj (|27|) is the 
thermal vacuum in the thermal BCS approximation. Hereafter, it will be denoted 
by |0(T);qp) CT . 

The average number of thermally excited Bogoliubov quasiparticlcs with quan- 
tum numbers jm in the BCS thermal vacuum |0(T); qp)cr is 

CT <0(T);qp|a] m a jm |0(T);qp) CT = y), (28) 

whereas the average number of nucleons in the same state is 

nj(T) = CT (0(T);qp|a 3 t m a im |0(T);qp) CT = u)y] + v)x) . (29) 

The function rij (T) determines smearing of the Fermi surface due to thermal and 
pairing effects. 

Equations (f2"2l with Vj , Xj (|2T|) are the well-known BCS-equations at finite tem- 
perature (see e^gpHUHEDl). 

Since the thermal vacuum |0(T);qp) CT contains a certain number of Bogoliubov 
quasiparticles, the excited states at finite temperature can be built on the top of 
|0(T);qp) ff by cither adding or eliminating a Bogoliubov quasiparticle. Due to the 
relations 

a] m |0(T);qp) a =a: i/ 5] m |0(T);q P ) CT , a^T); qp) r7 = a*y j ^ m .\Q(T); qp) CT (30) 

one can associate the first process with creation of a thermal quasiparticle having a 
positive energy, whereas the second process can be considered as creation of a tilde 
thermal quasiparticle having a negative energy. 

The simplest excitations on the top of the BCS thermal vacuum in an even-even 
hot nucleus involve two thermal quasiparticles. Their wave functions and energies 
are 

= J+) • 



[44 


>(r);qp)a, 


u) = e h +e 


[44 


>(r);qp) ff , 


w — e hh 


[44 


>(T);q P ) ff , 


u) = e h -e 


[44 




,,- 



(31) 

The square brackets [ in pip mean the coupling of single-particle momenta j\, 
j2 to the total angular momentum A with the magnetic quantum number fj,. 

As it should be, any thermal two-quasiparticle state with positive energy lu has 
a counterpart - a tilde-conjugated state with negative energy — u>. 

Quite interesting relations exist between electromagnetic transition probabilities 
to thermal two-quasiparticle states and their tilde counterparts. Let us write the 
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_EA-transition operator M. (EXfi) in terms of thermal quasiparticles 



-■ } x 



+ u< 3tl fay* [$Al£ - a *y^ ts^B-] 



(32) 



In ([31]), r^] a is the reduced single-particle matrix element of the .EA-transition 

operator; u£v = u n v n + u j2 v jl , i>Vj 2 = Uj x Uj 2 - Vj x Vj 2 and A = \J2X + 1. 

Only the terms containing the operators (jjj'a) and ^Xu(.?ij2) contribute to 
transitions from the BCS thermal vacuum state to any thermal two-quasiparticle 
state. Introducing the functions 

-1/2 



exp ( - ) - 1 



; X(u) = [1 + Y 2 (lo)] 1 / 2 



and taking the advantage of the relations 



(33) 



(34) 



we get the squared reduced matrix elements $| of the operator M. (EX/j,) between 
the BCS thermal vacuum and different thermal two-quasiparticle states (1311) 



iW .,(+) 



2^(+) 



2^(+) 



(35) 



From (|35|) it follows that transition probabilities to the states tilde-conjugated to 
each other and, correspondingly, having the energies ±cj relate with the factor 



*2M= exp (!)*£(-(«;). 



(36) 



Formally, the function Y(lu) in (|35[) is the Bose-Einstein distribution function which 
determines the average number of bosons with energy ui in a system in the thermal 
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equilibrium at temperature T . This gives an idea to treat two-fermion excitations 
in a hot nucleus as bosons. The probability to create a boson is proportional to 
the factor [1 + Y" 2 (w)], whereas the probability to annihilate it is proportional to 

Expressions (|35[) determine EA-strength distribution in a hot nucleus within 
the independent BCS quasiparticle approximation. An interesting point is that in 
contrast with T — case a portion of the -EA-strength appears in the negative 
energy region u> < 0, i.e., below the thermal vacuum state, at finite temperatures. 
This strength determines a probability of 7-ray emission by a hot nucleus, whereas 
the strength at u> > determines a photoabsorption cross section. Both the parts 
of the EA-strength contribute to the energy weighted sum rule (EWSR) at T 7^ 



EWSR = £ £ [*\([P]At) ~ 



T ii >h 



EE 

T ji >h 



+EE *u[*M^)-*wyi]i) 

T ji>h 

tl) 2 [*5£(«ffi) a o - vl - yl) - 4:lM;lm - 4)] ■ (37) 



3.3. Thermal phonons 

The second step in the diagonalization of the thermal Hamiltonian is to take into 
account the long-range particle- hole interaction i? p h- This interaction is respon- 
sible for the existence of different types of collective vibrations in nuclei. In this 
Subsection, the thermal quasiparticle random phase approximation in treating the 
vibrations of hot nuclei is discussed. 

Transformations (fT9|) and (|20|) with the coefficients determined by Eqs. pTj) and 
(|27p should be applied to the rest of the thermal Hamiltonian _ff p h — iJ p h as well. 
Then the thermal Hamiltonian (TT5|) takes the form 

^QPM = E E £ i(P)mPjm ~ Pjmftjm) 
t jrn 

~^E E (4 A) +p4 A ^^(r)M v (pr)-M^ (T) M v(/3r) }. (38) 

A^i rp=±l 

In terms of thermal quasiparticles the multipole operator Mj[ (r) has the same 
shape as the EA-transition operator ([32| . The only difference is that one should sub- 
stitute matrix elements r^j 2 for matrix elements /-^ 2 = (ji\\i x R\(r)Y\(6 , cf>)\\ j 2 ). 

The thermal Hamiltonian (|38p can be approximately reduced to the Hamiltonian 
of noninteracting bosonic excitations - thermal phonons. This occurs if one omits 
in (|38| the terms containing the operator B\^{jij2). Hereafter, the remaining part 
of the thermal Hamiltonian ([38| is denoted Htrpa- 

To diagonalize Htrpa, it is natural to use a trial wave function which is a linear 
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superposition of different types of thermal two-quasiparticle operators, namely 



T 3\32 



(39) 



The factors a and a* at the cross-ov er (i.e., ti lde-nontilde) terms of (f39|) are absent 
in the thermal phonon definition i n IIIHI2IIE] They appear d ue to ado ption of the 
new DTCR (JTUJ). The definition coincides with the one ir l 10 l 12 l 18 l when a = 1. 

One more very important assumption has to be accepted. We assume that the 
thermal biquasiparticle operators contained in (|39p commute like bosonic operators 



33^34'^-' 



b\\'&W'&3i33&323 



(40) 



and all other commutators are equal to zero. This assumption is known as the 
quasiboson approximation. 

Moreover, taking into consideration the long-range interaction we should rede- 
fine the thermal vacuum state. At first, we define it as the vacuum state |0(Q, Q)) a 
for thermal phonons 



(41) 



A thermal one-phonon state is constructed by acting on the thermal phonon vacuum 
|0(Q, Q))tr by the thermal phonon creation operator 



A thermal tilde one-phonon state should be define as 



yo(Q,Q)) a ee (-1) A -*QT_,J0(Q,Q)> C 



(42) 



(43) 



because just the operator (— 1) A ~ tl Q\_ fJji transforms under spatial rotations like a 
spherical tensor of rank A. 

The set of thermal one-phonon wave functions has to be orthonormalized. This 
demand together with the quasiboson approximation for commutators of biquasi- 
particle operators (|40|) imposes the following constraints on the phonon amplitudes 



1 T / 

\ " \ " a Xi -w Xl - +a Xi .w Xi . +f Xi - s Xi - ■ i x ' v Xi = fi S ; 

2 Z^Z^ yjl32 w 3l32 ^ yjl32 w 3l32 ^ L 3l32 a 3l32 T 1 '3l32 i 'jl32 ~ u " ' 



T 3132 



S^*S^ T a Xi - w Xl '- + a Xi w Xi ' +t Xi s Xi '- +t Xi - s Xl '- =0 

/ / / / y3l32 w 3l32 ^ y3l32 w 3l32 ^ L 3l32*3l32 x ''JUa °J1J3 U ' 



(44) 



T J1J2 
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In (H3J) , the following notation is introduced for the sums and differences of original 
phonon amplitudes: 

(l) =& ± &> (~) =fe ± &> 

^ ' hh ^ ' jija 

=rt h : (?) = ^ 2 ±^V (45) 

W ilia V*/ J U 2 

Moreover, ^ 2 ee (_!)*-*+* ^,^ a = (-1)*"*+* 

Constraints (|44[) imply that thermal phonon operators commute like bosons. 
With constraints (j44| one can find the inverse transformation to (|39| 

= E^hhQli + $L<h»i + + > 

i 

I = ^E^i^L + &» C Ki + VnnQU + %hQw> ( 46 ) 

i 

and then evaluate the following expression of the thermal RPA Hamiltonian Htrpa 
in terms of thermal phonon operators: 

HtRPA = ^2 E e 3 (PjmPjm ~ PjrnPjm) 
t jm 

~ I E JS E (*$ + ^i){[^<- 5 r^^'](Qi Mi +Q^)(Qt_ i; +Q A ^) 

Apii' rp— ±1 

+ [^5^-5^](Q^ i+ Q^)(Q^ /+ Q^,) _ (t.c.)} . (47) 

The notation " (t.c.)" in l|47p stands for the items which are tilde conjugated to the 
displayed ones. The functions D** and D^ 1 (r = n,p) are the following combinations 
of phonon amplitudes: 

ilia 

5 ? - " *VhVh&H) ^H7ly^A] ■ (48) 

ilia 

To find eigenvalues of WtrpAj we apply the variational principle, i.e., we min- 
imize the expectation value of Wtrpa over the thermal one-phonon state under 
constraints (|4~4"|) . 

It should be stressed that the phonon vacuum \0(Q,Q))<r is not the thermal 
vacuum in the sense of Eq. ^ and, thus, the expectation value of any physical 
operator with respect to \0(Q,Q))a does not correspond to the average over the 
grand canonical ensemble. 

After a variation procedure with respect to functions g, w, g, w, t, and s one 
gets a homogeneous system of linear equations. Since we use the separable effective 
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interaction, the system of equations looks simple 



1> 



1> 



Xi 



TJ1J2 

Xi 



TJU2 

Xi 



1 /* 



(A) ( + ) 



'0132 



E 

P =±i 



. .., ,_, /. ( K ) [ x 3i x h D 



Ai 
pr 



a 2 f {X) u {+) 



)Ai 



^i%2^pr] 



(7 XfaXjzDp,.] 



r 2 

J 3l32 3l32 



2A 2 ^T^f^ 



^ (4 A) 



y hXj2 D 



(49) 



'3132 1 ~' vt p=±i 

For a further discussion it is more convenient to rewrite the system (|4l?|) regard 
ing the functions and Z?p* as unknown variables. Then we get 

Xi / X 

1 —l / D 



5> 

p=±l 



(A) 




= [^w] _1 



D 



where the function X^(u>) reads 

1 V- T / f (A)^2 



A 2 



J1J2 



(f- ■ r 



(»ai) 2 ^d-^- 



y 2 ) fw ( " ) ) 2 eH (V -y 2 ) 



(+K 2 



^1.72' 



(50) 



(51) 



Demanding the existence of a nontrivial solution to the system of linear equa- 
tions (|50p , we derive the secular equation for the energy of the thermal one-phonon 
state u>\i 



1. 



(52) 



Equation ([52]) is strictly the same as was obtained by other methods j n j 31 | 3 2 | 3 3 j| . 
does not differ as well from the previous results of the TFD-QPM approach! ^^ 

One obtains the same secular equation as (I52p if the variational procedure is 
applied to the expectation value of Htrpa over the tilde thermal one-phonon state 
Ql_.|0(Q, Q))a- As it was stated in the previous subsection, both the positive and 
negative energy excitations of a hot nucleus built on the top of the thermal vacuum 
have a physical meaning and should be considered on equal footing. So the negative 
roots of (|52|) are identified with energies of tilde-phonons. 

Unfortunately, the values D$ l and D^ 1 cannot be determined unambiguously 
from equations (f50|) and the normalization condition ([44]) . The reason for this in- 
determinacy is the following. Equations ([5H|) enable one to prove that the TRPA 
Hamiltonian 7Ytrpa is diagonal in terms of the thermal phonon operators, i.e., 



HTRPA = y2^Xi(QxaiQxai ~ QxuiQxui) 



Xfli 



(53) 



It is seen that the Hamiltonian (|53[) is invariant under the unitary transformation 
which mixes non-tilde and tilde thermal phonons but preserves the secular equation 
(f52|) . Just due to this invariance the systems of equations (|49|) or (|50|) cannot be 
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solved unambiguously. To overcome this problem, one needs to involve additional 
considerations. To this aim, minimization of the thermodynamic potential was used 
while considering the pairing correlations at finite temperature (see Subsec.( |3.3l) ). 
However, in the case of a thermal phonon system this procedure is not so straight- 
forward. 

Let us turn back to equations (|49[) . The following interrelations between D^ 1 
and D^ 1 can be easily found from the normalization condition (|4"4l : 



D 



Xi\ 



D 



(54) 



where A/" Al is given by 



A'' 



Xi 



d_ 



'i-x£ A) ( WAi )(4 A) +^ A Y 



(A) 



(55) 



Instead of the functions D^ 1 and D^ 1 it is convenient to introduce the new 
variables 



2A 2 X< A) (w A< ) 



1a 



2X 2 X^(u Xi ) 



(56) 



which obey the following condition: X^ — = 1. 

Then we substitute X\i and Y\i for D^ 1 and D$ l in equations (|49|) and get 



1> 



Xi 

3132 

Xi 

3132 



1 



f (A) U ( . + ? 
■>3l32 JU2 



(x h x h X X i + o- 2 v n y j2 Y Xi ) 



f (A) u {+) 

J 3132^3132 I y 

\V3lVh X Xi 



(7 X O-. X 'icy 



j 1 Xj 2 I\i) , 



A; _o „(-) 



j;(A) (-) 

Jjl32 V 3l32 ( V . ,, ,. \ . \ 



(57) 



Some conclusions concerning variables X\i and Y\i can be achieved if one an- 
alyzes their behaviour at the limit T — > 0. Obviously, the phonon wave function 
([39f at T — should consist of two components only: a'^ aj 2 and aj 2 aj 1 . Since at 
T — > the thermal occupation numbers of quasiparticles yj tend to zero (xj — > 1), 
the demand can be fulfilled only if y A i — > (and synchronously X\i — > 1). In such 
a case the tilde-amplitudes and A jj 2 tend to zero whereas ipj*j 2 and 4^1 j 2 

survive. 

Let us define the thermal phonons with the amplitudes corresponding to Y\i — 
and X\i — 1 as gL, g^^i and name them "reference phonons" (r-phonons). The 
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vacuum for r- phonons is denoted by \0(q,q)) a . From (|5T[) one can conclude that 
the thermal phonons corresponding to other values of X\i , Y\i can be produced by 
applying the unitary {X\i, Y\i} transformation to the r-phonons 

QU = ^U~Y^- (58) 
The vacuum for the phonons Q! j, Q| . has the forrrJU 

|0(Q,Q)) CT =exp{- J ftr 6 /2}exp{^ 9 i^ i }|0(g,g)) (T , (59) 

where Kb is the entropy operator for noninteracting bosons 

Kb = ~ E<«U«A* lnY ^ lnX AJ- (60) 

Now we are ready to determine the values X\i, Y\i- To this aim, we minimize 
the thermodynamic potential Qb for the r-phonon system which reads 

n 6 = Ao(Q,Q)\J2 u, ^U^i ~ TK b \o(Q,Q)) a 

Xfii 

= {ujmYI + T(Y* In Y* + Xl \nXl)} . (61) 

Xfti 

Varying (|6 1 (1 with respect to Y\i and equating the result to zero we get 

Y Xl = [exp (^) - l] _1/2 , X Xt = [1 + Y^ 1/2 . (62) 

The coefficients Y Xi are the thermal occupation factors of the Bose-Einstein statis- 
tics. They determine the average number of r-phonons in the thermal phonon vac- 
uum 

AQ(Q,Q)\ QlnQxtfl o(Q,Q))* = Yxi • (63) 
It is worthwhile to note that the problems of correct construction of the thermal 
RPA phonon operator and the thermal RPA vacuum state were already discussed 
infusing the BCS Hamiltonian as an example. The authors ol^also confronted with 
an ambiguity of RPA solutions at finite temperature. To get it out, they introduced 
the " non- rotated" thermal pairing phonons and then made a thermal rotation which 
minimized the thermodynamic potential. At the same time, these circumstances, 
i.e ., the invariance of the Htrpa under the thermal rotation (|58|) . were overlooked 
i n | 10 | 12 | j n t nese papers, it was assumed that Y\i = at any temperature and the 
vacuum for "r-phonons" was identified as the "true" thermal vacuum state. That 
is why ther mal boso nic occupation numbers did not appear in the corresponding 
formulae in! 10 ' 12 ! 18 ! 

Now we can fix the factor a. To this aim, we consider the behaviour of thermal 
phonon amplitudes (157p when the coupling constants of a separable multipolc in- 
teraction Kq^I tend to zero. When n^l — > the thermal phonon energy uj\i — ► £^j 2 
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and only one amplitude has to survive in the corresponding phonon wave func- 
tion. In particular, all the amplitudes of backward going components vanish when 
Kq A i ~^ 0. Taking advantage of (|34| it can be shown that if uj\i — » e^. ^ £,^j 2 ~ * 
whereas (Vjij 2 ) 2 — > 1 at any a value. 

When uj\i — > £^- a the amplitudes <j> and tp also vanish at any value of a. However, 
limiting values of tp and <fi appear to be dependent on a, namely, 



lim _ x n X J2 X (s^l ) + a 2 y h y j2 Y (e)^ ) _ y2/ j+) ^ „ 2v 2 f J+) 



_(+) 



lim 



^-vl-vl) 1 ' 2 



X^y^^ + l). (64) 



Since from a physical point of view 4>jlj 2 should vanish, we choose a = i. Then not 
only 4>f ih = but also ipg ja = 1. The present result on DTCR agrees with the 
conclusion inP^ Inj^ the appropriate choice of DTCR was specified by the fermion 
number conservation in the system. If the number of fermions in the system is not 

conserved, DTCR should have the form ([9]), i.e. A = A. This seems to be just our 
case since the number of quasiparticles in a nucleus is not conserved. 

Now we complete constructing the thermal phonon operator and the "true" 
thermal vacuum state in the random phase approximation. They should be calcu- 
lated with formulae (f39|) and ([57)) where X\i and Y\i are given by ([62]) and a = i. 
Hereafter, the thermal phonon vacuum is denoted by |0(T); TRPA). We would like 
to stress once again that the thermal RPA vacuum |0(T);TRPA) reduces to the 
thermal BCS vacuum when the particle-hole interaction vanishes only if the coef- 
ficients X\i, Y\i are the phonon thermal occupation numbers given by (|62p and 
a = i. 

At the end of this subsection we calculate a matrix element of the EX- transition 
between the RPA thermal vacuum and a thermal one-phonon state. To this aim, 
one has to write the operator Ai(EXfi) (I32| in terms of thermal phonon opera- 
tors. Taking into account only the term of ([52")) with the operators A\^(jij 2 ) and 
A \jl(hh) one gets 



% T j ± j 2 

x { [ u hU x n x j*$ h + VhVh9^ j2 ) + ^H^hvA] (QU + QxHi) 

= I E E T r {Xm (Q[, t + Q-J + Y Xl (Qt_. + Q }, (65) 
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where 



71 la V-' V T 



Cu (+) ) 2 e (+) fl- v 2 - v 2 ) (v y ~') 2 e { ~> (v 2 - v 2 ) 



(66) 



Expression ([63)) differs significantly from that mP^Sl Specifically, there appear 
terms with tilde-phonon operators in ([65]) . This is a consequence of {X, Y } rotation 
of thermal phonons. The item proportional to the factor Y\i is responsible for 
transitions to tildc-phonon states lying below the thermal vacuum state, i.e., for 
decay of the thermal vacuum. 

Thus, there are two types of matrix elements corresponding to excitation and 
de-excitation processes of the thermal vacuum 

<f Aj = {Q{T) ] TR¥^M{EX l i)Q{j{){T) ] TRVK)=X x ^Y^\ 

T 

$ Aj = (0(T);TRPA||M(^A/*)Q^.||0(T);TRPA) - Y Xl ^r- (67) 

r 

The factors X\i and Y\i in (|6"T)) were missed iJI^HBI since the r-phonons q\ ■ , q\^ 

were explored rather than the "rotated" Q\ ai ,Q\ ai - 

The factor X 2 t in the photoabsorption cross section for hot nuclei as well as 



1041 

its role were discussed in detail in!^ This factor occurred also in the response 



function i djfj bu t was missed in the B(EX) j expression inH^ Note that in both 
the papers^2123| the thermal RPA equations were derived with the equation of 
motion method for bifermion operators a] a\ and a\ a.,-, . 

1 Jl ]2 01 J 2 

As in the case of transitions to thermal two-quasiparticle states (see Sub- 
sec. (1321)), there exists the following relation between transition probabilities from 
the non-tilde and tilde one-phonon states 

$ 2 Aj - exp(aWT)$L (68) 

This relation is equivalent to the principle of detailed balancing connecting the 
probabilities for the probe to transfer ene rgy u> to the heated system and to absorb 
energy to from the heated system (see e.gH^). 

The model energy-weighed sum rule at the thermal RP approximation is given 

by 

EWSR = ]T WA i(*^ - $L) = E w * ( T P + T nf ■ ( 69 ) 



It is worth-while to note that the EWSR (|69|) appears to be independent on the 
thermal phonon occupation factors. The numerical value of (|69| should coincide 
with that of (l37| . 

[00] 

Nominally, this expression coincides with that in. 1 ^ However, the essential dif- 
ference between the present result ancP^ is the contribution of TRPA states with 
negative energies ui\i (i.e., tilde states). Inp^the total strength of _EA-transitions is 
located in the positive energy region whereas in the present approach some fraction 
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of strength is carried by the states with ui\i < 0. As a result, part of the EWSR per- 
taining to positive excitation energies is greater than the total one that is obvious 
from (1691). 



3.4. Interaction of thermal phonons 

In this Subsection, we deal with a coupling of thermal phonons, i.e., go beyond the 
thermal RPA. Physical effects which can be treated in this order of approximation 
relate to fragmentation of basic nuclear excitations like quasiparticles and phonons, 
their spreading widths and/or more consistent description of transition strength 
distributions over a nuclear spectrum. The problem of temperature dependence of 
the gia nt resona nce width which was intensively discussed not long ago (see, e.g., 
reviewJS^SSEZlj a i so belongs to this set. 

The part of Ti. (fT8"|) which is responsible for these effects is the so-called 
quasiparticle-phonon interaction (or the cubic anharmonic term) 7i qp h 

1 f (A) 
= "a E EE T ^f + Qx#)***(*h) + (h-c.) - (t.c)} , 

Xfli T j ± j 2 V T 

-4:1 (*& i4/y i + yfh, $kh£) ■ ( 7 °) 

The coefficients Xj*j , D^j 2 and Zjfj 2 are given by 

*\ Xl fx\ fy\ m , s 

yj = x jl x j2 y y j + y n y l2 y x j , Z jlh = x h y h X\i + y n x l2 Y Xl . (71) 

The interaction 7i qp h couples the multiphonon states whose structures differ by 
one phonon, i.e., one-phonon states with two-phonon states, two-phonon states 
with three-phonon states etc. This term can be named the "cubic anharmonic" one 
because in the lowest order of boson expansion the fcrmionic operator i?A^i(jij2) is 
substituted by the product of two bosonic operators (or phonons). The remaining 
part of the thermal Hamiltonian consists of the items ~ &X— uO'iia)^*, u 0*3.74) which 
are equivalent to the sum of products of four bosonic operators. Its contribution 
will not be considered here. 

To take into account the term 7i qp h, we again apply the variational principle. 
A new trial wave function is assumed to be of the form 

I M 



\9„(jm)) = {E^( J ^M i +E p A 2 t(^)[QL 1 «LJ: 

]S X x± (J«0 [QUQ'-J m + J2 T ^ W [QhiPhJ i }|0(T); BPA), (72) 
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where R, P, S, T are the variational parameters which should be determined. As 
one can see in ([72]) . the thermal vacuum is kept the same as in the TRPA. At T = 
the latter approximation is valid if the quasiparticle-phonon interaction is relatively 
weak. 

The trial wave function (|72p has to be normalized and, therefore, the amplitudes 
R, P, S, T should satisfy the following constraints: 

£[W)] 2 + £ {2[P^(Jvtf + (■/")] 2 + 2[T^(^)] 2 } = !• (73) 

i Aiii 

^2*2 

Since the trial wave function contains three different types of components, there 
are three types of interaction matrix elements which couple a thermal one-phonon 
state with two-phonon ones 

( J*) = ( RPA ; °( T )\Q JMiH qph [Q{ in Q{ 212 ] J M \0(T); RPA) , 
V££(Ji) = (RPA;0(T)|Q JMi W qph [Qt iii Qt_.ji f |0(T);RPA), 
W££(Ji) = (RPA;0(T)|Q JMi W qph [Qt_. i Q+-.J^|0(T);RPA). (74) 

The matrix element U^ 1 ^ (Ji) correspond to a transition from one to two phonons, 
whereas the matrix elements V^*^ (Ji) and W^^fJi) describe a thermal phonon 
scattering and a thermal phonon absorption, respectively. The explicit forms of 
U, V and W are given in Appendix A. 

Thus, we should minimize the expectation value of the thermal Hamiltonian 
Wrpa +T~Cq P h over \^ V (JM)) at constraint (|75|) . The expectation value is given by 



(tf„(JAf)|WHPA + Wqph|*«(JM)) =^2uji[Ri{Ju)] 

i 

+ 2 h,, + "x^)[Px£ w] 2 + E - ^ ) K^ 1 ( J ^)] 2 

Xlil i Xiii 

A2«2 ^2«2 

+ S^(Ju)W^(Ji) + T^(Ju)V^(Ji)}. (75) 



After the standard variation procedure one gets a homogeneous system of linear 
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equations (fy, is the energy of the state \$> V (J M))) 



+ S^(Ju)V^{Ji)+T^{Ju)W^{Ji)} 
(Ji/)(« Al<1 + w A2i2 - trf + i £ Ri(Jv)U&£ (Ji) 



o , 







. 



. 



(76) 



The system has a solution if r\ v is the root of the following secular equation: 



Minimizing the expectation value of the thermal Hamiltonian 7iRPA+Wqph over the 
wave function \^f v (JM)) we get the corresponding equations for tilde-conjugated 
states. They also can be obtained from ([76)1 and (|77|) by changing the sign of r\ v . 
That is if r\ v is the energy of a given state |$ „(J M)) , then — r\ v is the energy of the 
tilde-conjugated state \^ V {JM)). 

Certainly, the above results deviate from those of ff^ Due to the changes of ther- 
mal RPA-phonon amplitudes and, in particular, their dependence on the thermal 
phonon occupation numbers, the quasiparticle-phonon interaction (|70|) couples one- 
phonon states with three types of thermal two-phonon states whereas i only one 
type was taken into account (namely, the first term in (|72|) . 

The terms Q^Q^ and Q^Q^ in the trial wave function describe the processes 



which were not considered irM^ and later inM^ The processes of thermal phonon 
scattering and absorption became possi ble d ue to the presence of thermal r-phonons 
in the thermal phonon vacuum. Since irP^ the thermal quasiparticle-phonon inter- 
action was treated on the basis of thermal vacuum for r-phonons, the above effects 
were missed. 

The inclusion of the terms Q^Qt and Q^Qt in the trial wave function produ ces 
the new poles in the secular equation (|77p in comparison with the previous studyf^ 
namely, (ojx^ -wa 2 . 2 ) and -(w Alil +u\ 2 i 2 ). 

The new equations of the TFD-QPM approach (|76l) and (|77|) are in qualitative 




(77) 




agreement with those i and. i 



in some respects, in. 1 
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As it was mentioned in the Introduction, the consideration h i I I was based on 
the Matsubara Green function technique and Nuclear Field Theory. The approach 
oP21 i s especially close to our approach since it treats a hot nucleus as a system of 
interacting TRPA phonons. Although the equations of both the approaches seem to 
be hardly compared "term-to-term" because of quite different formalisms explored, 
in both the cases the negative TRPA roots are in the game, the poles in the equa- 
tions are of the same types, phonon interaction matrix elements are similar etc. 
That is why we can establish their "qualitative" agreement. 

The formulae for the matrix elements of EX transitions from the thermal vacuum 
state to the state (|72|) and its tilde-counterpart are quite obvious. In the leading 
approximation the operator M. (EX) induces a transition from the thermal vacuum 
to onc-phonon components of the thermal state Y$> V (JM)) (or |\P „(J M))) only. 
Matrix elements of direct transitions from |0(T);RPA) to two-phonon components 
are very weak although do not vanish. Thus, one gets 

$(JV) = (^ u (JM)\\M(EX)\\0(T);RPA) = ^^(Jz/)^, 

i 

i(Ji/) = (^(JM)||M(£A)||0(T);RPA) = J^M^Ji- (78) 

i 

where $ji and <E>ji are given by (|67|) . 



4. Summary and conclusions 

The pre sent stud y was motivated by the necessity to reexamine the TFD-QPM 

a p prc .adp™Hr in theory of hot nuclei. In this respect, a couple of new facets 

of the TRPA formulated within the TFD were found and their influence on the 
coupling of thermal phonons was established. 

We showed that for the Hamiltonian consisting of a mean field, the BCS pair- 
ing interaction and separable particle-hole effective interactions, amplitudes of a 
thermal phonon wave function (|39p were not determined unambiguously by diago- 
nalization of the RPA-part of the thermal Hamiltonian. To fix the coefficients of a 
linear transformation from the set of thermal two-quasiparticle operators to phonon 
operators, one should impose an additional demand — to minimize the thermody- 
namic potential of the system of free thermal phonons. This was achieved by using 
one more unitary transformation — a thermal rotation of " reference" phonons (r- 
phonons). As a consequence, the Bose-Einstein thermal occupation numbers (ther- 
mal occupation numbers of phonons) come to play. It should be stressed that the 
new ingredient does not affect the main TRPA equation, i.e., the equation for ther- 
mal phonon energies. 

As far as we know, the thermal rotation of phonons in the TRPA framework was 
discussed only irP^ when treating the pairing BCS Hamiltonian. The ther mal rota- 
tion of phonons was missed not only in the preceding TFD-QPM 

studies^™ 

but it was not also mentioned inP^ Moreover, the Bose-Einstein thermal occupa- 
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tion numbers do not appear when the TR PA equa tions are derived by exploring 
the phonon operator with scattering termsPSEHEl] 

Thus, according to the present results amplitudes in the thermal phonon wave 
function depend not only on thermal occupation numbers of Bogoliubov quasi- 
particles making up the phonon but also on the thermal occupation numbers of 
the phonon. Moreover, presently the corresponding thermal phonon vacuum con- 
tains some amount of r-phonons with probabilities determined by their energies in 
accordance with the Bose statistics. Due to this, the processes of excitation and 
deexcitation of a hot nucleus can be regarded on equal footing. 

The above results have weighty consequences. The first one is the appearance of 
the factor 1/(1 — exp(— w/T)) in the EX transition strength. Its role was discussed 
in34 

where it was derived using the Green function technique (see alscE21). This 
factor somewhat enhances the low energy part of the EX transition strength. It is 
also essential when to can take negative values, i.e. when considering a decay of a 
hot nucleus. 

The second consequence concerns the thermal quasiparticle-phonon interaction 
Wqph- This consequence is two- fold: a) renormalization of phonon-phonon inter- 
action vertices (cf. matrix element U^^(Ji) (jA.lj) with that irP2); b) essential 
extension of the thermal two-phonon configuration space and the corresponding 
complication of trial wave function (|72l . 

With the above new ingredients the TFD-QPM approach well conform (at least 
qualitatively) with a more traditional one exploring the Matsubara Green function 

|01 |99| 

method! 1 'Now in both the approaches one has two kinds of thermal occupation 
numbers (fermionic and bosonic) and the identical sets of processes giving rise to 
the phonon-phonon coupling. 

In view of the above discussion, we conclude that the results of calculations 
irffSl should be revised. At the same time, it should be stressed that the thermal 
rotation of phonons cannot be regarded as a mandatory ingredient of a microscopic 
treatment of boson-like excitations in many fermion systems. Its necessity seems 
to b e intimately re lated with a two-body interaction character. For example, in 
RefsJEEIE^ESl some approximations going beyond the TRPA were constructed 
by using the TFD formalism and their validity was examined by the example of the 
Lipkin model. In those studies the transformation from thermal bifcrmion operators 
to thermal phonon operators was unambiguously determined by diagonalization 
of the thermal model Hamiltonian. No additional demands or assumptions were 
needed to be involved. 

One more new feature of the present study is exploring the double tilde- 
conjugation rule in the form proposed by I. OjimaPSl More precisely, we found 
that just the Ojima's formulation of DTCR guarantees the correct behaviour of the 
thermal phonon wave function in the limit of vanishing particle-hole interaction. 
Furthermore, the Ojima form of DTCR implies the complex thermal rotation for 
fcrmions. An interesting point is that the effect of the DTCR choice was revealed 
only at the TRPA stage while constructing a thermal phonon operator, whereas at 
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the pure fermionic stage when, e.g., pairing correlations were considered, the role 
of DTCR did not show up in and both the versions gave the same results. 
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Appendix A. 

Here we display the matrix elements of the quasiparticle-phonon interaction 7i qp h 
between thermal one-phonon and different two-phonon configurations. The phonon 
amplitudes are shown in (|57p and include both quasiparticle and phonon occupation 
factors. The symbol < ■ ■ ■ > is the standard 6j-symbol. 



n-Ji (■ ; ;\ yJi „,( — ) f./.^ih A X 2i2 i jAiii iA 2 i 2 . Aiii t A 2 i 2 i fAi ii ~A 2 i 2 \ 

- y J \ W H U Xlll 4> X212 + Z Xlll ri X2i2 +« Alil f A2i2 +f Xili ri X212 ) 

Ml 3l32\^ 3233 ^ 3331 V 32]3 V 3331 '3233 ^3331 S J2J3 '3331 I 

+ Z' H ?/ + ' (J, Xlil ? X2i2 + J A i*i^ A 2i2 _ iA 2 i 2 _ ?Alil„/.A 2 i 2 \ 

3132 jlh \ V 3233 ^3331 v 3233 '33 jl '3233 ^ 3331 ^3233 ^ 3331 I 

- Z Ji ?/ + ' U Alil f A2 * ! -\- (h Xlil n X2 ' i2 — n Xlil Z X ' 2i ' 2 — C X i i iJ, X 2i2\ 

3231 3132 \V 3233 ^3331 ^ v 3233 <333l 'j'aj's ^j'sj'i ^3233^3331 > 

- y Xlll v { ~ ] U X2l2 *b Jl - + 6 X2l2 (b Jl - +rj X2Z2 n Jl - +£ X2i2 £ Ji .) 

^3132 3132 \ ^3233 v 333l ~ ^ 3233 v 333l ~ '3233 '3331 ~ ^3233 ^3331 1 

+ 2 Al ")/' + ' fi? 2 ' 2 ^' +7) X2i2 f Jl —n X2i2 tb Ji —J X2i2 d> Ji \ 

3132 3132 \ y 3233 >3331 T v 3233 ^3331 '3233 * 33 jl ^3233 v 3331 1 






(A.l) 



where 




September 8, 2008 8:59 WSPC/INSTRUCTION FILE dzhiocv 



24 Authors' Names 



V, 



Aiii 



(ji) = E 



3U233 
c(Ai) 



f (-0 
■'3132 



J Ai A 2 
J3 ji h 



+ 



_f_]\Ai+A 2 + J 



Ai A 2 J 
33 ji 32 



Kxl? 2 ji(3ih33) 



(A 2 ) 



A2^2 



A 2 Ai J 
j3 ii h 



F\rZji(3i3233) 



(A.2) 



where 



cJi I a a a \ yJi ,n( — ) f»/iAi*i A x 2i2 I j.Ai»i / A2»2 I „Aiii t A2i 2 _i_ fAin A2i2\ 

- y Jl - (ib Xlil d> X2i2 + 4> Xltl tb X212 +n Xlil f X2i2 + f x ^n x ? i A 

3\32 u 3l32 \^3233 Y 333l ^ v 3233 ^333l '3233 ^3331 ^3233 '3331 > 

I ■prJi ,.(+) (7Miic^2i2 I JAiii A 2 »2 _~Aiii7A 2 i2 _ ?Aiii TA 2 i 2 \ 

^ ^'U2 "jlj 2 V^2i3 ?J3J1 ^ ^2J3 'O'ajl ''.72.73 ^jsj'l S'2J3 ^J'SJI ) 

— Z" Ji ?/ + ' ) (ih Xlil 7 X2i2 JL f / ) Aiu~A 2 i2 _„Aiii /A 2 i 2 _ tAiii / A 2 i 2 \ 

■72.71 3l32\ v 3233 ^3331 ^3233 '3331 ' 3233 V 3331 ^3233 ^3331 / 



ri&jtiM - -"jIMj, -^jtnL 

- V Al4l w ( ~ ) f-(/> A2l2 ^ j4 - + 6 X2l2 d> H - +n X2l2 n Jl - + f X2l2 f Jl . ) 

^3132 3132 \ ^3233 ^3331 ~ v 3233 v 333l ~ '3233 13331 ' ^3233 "=73.71/ 
1 -p-Aiii (+) / / A 2 i 2 ~Ji _i_J,A2i2?/i _ „A 2 i 2 / Ji _ t^2i2 xJi \ 

^ JU2 "jl.72 1^2.73 'Osjl ^ ^j 2 j 3 ?J3J1 '.72.73 V 333l ^3233 ^333l) 

_ ^Aiii (+) /7A 2 i 2 Ji _i_jA 2 i2fJi _ JfA 2 i 2 Tji _7A 2 i 2 X/i \ 

.72 71 ".71.72 V^j 2 .73 'J3J1 ^ ^j 2 j 3 ^J3J1 '.72.73 ^73 jl S7273 ^3j'l / 



.Fr 2 * 2 , f 7i 797,) = X X2i2 v^ (ib Xlil tb Ji - +d Xlil 6 Ji - +ri Xlil n Ji +£ Xlil f Ji .) 

■'Aiii Ji\JU2J3) ' x j 1 j 2 u 3l32 \r3233 V 333l ^ V 3233 ^3331 ^ '3233 ''3331 T S273 ^333l) 

- y X2i * v (~} (ih Xlil th Ji . + 6 Xlll (b M - +n X2i3 n Ji - + f X2i3 f Jl ) 

^3132 3132 \ ^3233 ^3331 ~ v 3233 v 333l ~ '3331 '3331 ' ^3331 ^3331 ) 
1 ^?A 2 i 2 (+) / , Aiii Ji 1 _ „Aiii Tjj _ \ 

' 3132 3132 VT3233 '3331 ' ^3233 ^jsjl '3233 ^3331 ^3233 "333l> 

_ p-A 2 i 2 ( + ) (JMii^Ji +7> Xlil 7 Ji —r) Xlil il> Ji — 7 Xlil th Ji ) 
3231 3132 \ v 3233 13331 ^ v 3233 ^3331 '3233 ^3331 ^3233 v 3331 1 ' 



313233 

f(Ai) 
3132 



33 ji 32 



f 



_^_^^Ai+A 2 +J 



Ai A 2 J 

33 jl 32 

f (A 2 ) 
J 3132 

A 2 l2 



n \2i 2 Ji(3i3233) 



A 2 Ai J 

33 jl 32 



n X2 X. H {3i3233) 



(A.3) 



where 
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Qi { x ■ ( ii iiii) = X Ji - v^~) (ih Xlil S X2i2 + S Xlil ib X2i2 +n Xlil £ X2i2 + £ Xlil n X2i2 ) 

- y M uH (ib Xlil 6 X2i2 + d> Xin ib X212 +n Xlll £ X212 +f Xin 7j X2i2 ) 
■S]i]2 3132^3233^3331 ^3233^3331 13233^3331 ^3233 ' 

i v ,n ?/ + - ) (ib Xlil f X2i2 +(h Xlil n X2i2 — n Xlil 6 X2i2 — £ A i i i ? p V2i2N 

— ? Ji ii^ +S> (7, x i i i'f X2i2 -L A> Xlil n X2i2 — ri Xlil rh X ' 2i ' 2 — 7 Xlil i// 2 * 2 N 

hn^hn^hh^hh ' ^j 2 j 3 'hh 'iih^hh ^hh^Hh' 



- V Xlil v^ U X2%2 ib Jl - + 6 X2l2 d> H - +n X2i2 n Ji - +£ X2t2 f Jl . ) 

3132 3132^3233 v 333l ~ V j2j3 ^3331 ~ 'I3233 '3331 ~ ^3233 ^3331 > 

_i_ -7*1*1 (+) (7\ 2 i2 Ji 1 7^2 fJi _7y X 2i2jJi _ ?A 2 i 2 7 Ji \ 
' 3132 3132 ^3233 13331 ~ * 3233 ^3331 1 33 jl ^3233 v 333l> 

— 2 Alil ?/' + ' (ib X2i2 r) Ji A- (h X2i2 T Ji —n X2i2 ib Ji —F X2i2 (h Ji ^ 

32 ji 3132W3233 '3331 ^ ^ 3233^3331 13233^3331 ^3233 V 33 ji 1 • 
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